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We investigate the baryon mass spectrum in quenched lattice QCD with exact chiral 
symmetry. For 100 gauge configurations generated with Wilson gauge action at /3 — 6.1 on 
the 20 3 x 40 lattice, we compute (point-to-point) quark propagators for 30 quark masses in 
the range 67 MeV < m q < 1790 MeV. For baryons only composed of strange and charm 
quarks, their masses are extracted directly from the time correlation functions, while for 
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those containing u(d) light quarks, their masses are obtained by chiral extrapolation to 
m n = 135 MeV. Our results of baryon masses are in good agreement with experimental 
values, except for the negative parity states of A and A c . Further, our results of charmed 
(including doubly-charmed and triply-charmed) baryons can serve as predictions of QCD. 

1. Introduction 

One of the basic objectives of lattice QCD is to compute hadron masses nonpertur- 
batively from the first principles. For hadrons only composed of strange and/or charm 
quarks (e.g., in Fig. [IJ, their masses (in quenched approximation) can be measured 
directly with presently accessible lattice sizes. However, for hadrons containing u, d light 
quarks, the performance of the present generation of computers is still inadequate for 
computing their masses at the physical scale (m n ~ 135 MeV), on a lattice with enough 
number of sites in each direction such that the discretization errors and the finite volume 
effects both are negligible comparing to the statistical ones. Nevertheless, even with lat- 
tices of moderate sizes, lattice QCD can determine the parameters of the hadron mass 
formulas in the (quenched) chiral perturbation theory. Then one can use these formulas 
to evaluate the hadron masses at the physical scale. In this paper, we extrapolate the 
baryon mass linearly in m\ (e.g., S+ in Fig. EJ). 

Here the quark fields are formulated in the framework of optimal domain-wall fermion 
[d such that the quark propagator and its fermion determinant can attain the maximal 
chiral symmetry for any finite N s . For 100 gauge configurations generated with Wilson 
gauge action at (3 — 6.1 on the 20 3 x 40 lattice, we compute (point-to-point) quark 
propagators for 30 quark masses in the range 0.03 < m q a < 0.8 [|2J- Then we determine 
the inverse lattice spacing a -1 = 2.237(76) GeV from the pion correlation function, with 
the experimental input of pion decay constant f n = 132 MeV. The strange quark bare 
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Figure 1. The even parity amplitude A+(t) 
of the time correlation function of is plot- 
ted versus the time slice. The solid line is 
the single exponential fit for 10 < t < 15 
where the effective mass attains a plateau, 
and it gives M n o {1/2+) = 2677(30) MeV, in 
good agreement with the mass of f2°(2697). 



Figure 2. The mass of the doubly charmed 
S+ baryon versus the pion mass square, 
for J p = l/2 ± states respectively. The 
solid lines are linear fits using the small- 
est 11 masses with 0.03 < m^a < 0.08. 
The even parity state gives Mg+pM+j = 
3522(16) MeV, in good agreement with 
H cc (3520) observed by SELEX [EJ. 



mass m s a = 0.08 and the charm quark bare mass m c a = 0.8 are fixed such that the 
masses extracted from the vector meson time-correlation function are in close agreement 
with the experimental masses of 0(1020) and J/?/>(3010) respectively. Then the masses of 
any other hadrons containing c, s,u and d quarks 2 are predictions of QCD from the first 
principles. 

2. Baryon interpolating operator and time-correlation fnnction 

First we define the notation for the "diquark" operator 

r A r B]_ F /AB B A \ 

n 1 CL \xa = Cabc-L a/3{<l xabQ x/3c ~ <l xabQ x/3c) 

where q A and q B denote quark fields of flavors A and B\ e a b c is the completely antisym- 
metric tensor; x, {a,b,c} and 7} denote the lattice site, color, and Dirac indices 
respectively. Here T a p = —Yp a such that the diquark transforms like a spin singlet. 
Further, we define 

(q A rq B ) x . a = £ a bcT a f3<i A xabCl 3 X f3c 

2 In this paper, we work in the isospin limit m u — rrid- 
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Then the interpolating operators for N and A can be written as 
N X1 = [u(C7 5 )d] xa u X7a , A X1 = (dC , 7 M d) a;a d :E7a 

where C is the charge conjugation satisfying Oy^C -1 = — 7J and (C*7 5 ) T = — C75. 

Now suppressing the Dirac and site indices, the interpolating operators for baryons 
with J = 1/2 can be written as 

N = [u(C 75 )d]u, 2° = [u(C 75 )s]s, £+ = [ U (C 75 )s]u 
A = [d(C 75 )s]u + [s(C 75 )u]d - 2[u(C 7 5)d]s, 

while the interpolating operators for baryons with J = 3/2 can be written as 

A~ = (dC 7M d)d, Q~ = (sC 7At s)s, £+ = (uC 7At s)u + (sC 7At u)u + (uC^ujs, 
2° = (uC7 M s)s + (sC , 7 Ai s)u + (sC7 M u)s 

Similarly, the operators for charmed baryons can be constructed, e.g., = [c(C7s)s]s, 
and 3+ = [c(C T5 )d]c. 

The time-correlation function of any baryon operator B is defined as C a p(t) = J2x(B xa Bop) 
which can be expressed in terms of point-to-point quark propagators from y = (0, 0) to 
x = (x,t). Its average over gauge configurations is fitted by the usual formula 

(C7(i)) = lipi(Z +e - m+at - z_e- m - a{T ~^) + l^J±(Z + e- m + a ( T -Q - Z_e" m - ai ) 

where m± are the masses of even and odd parity states. Thus, one can use parity projector 
(1 ± 74) /2 to project out two amplitudes, 

A+(t) = Z + e- m+at - Z_e" m - a(T -' } , A.(t) = Z + e- m + a{T ' t) - Z_e~ m - at . 

Now the problem is how to extract m± from A± respectively. Obviously, for sufficiently 
large T, there exists a range of t such that, in A±, the contributions due to the opposite 
parity state are negligible. Then m± can be extracted by a single exponential fit to A±, 
for the range of t in which the effective mass m e ff(t) = \n(A±(t) / A±(t + 1)) attains a 
plateau. On the other hand, if T is not so large, then it may turn out that the heavier 
mass, say m_ (assuming m_ > m+), could not be easily extracted from A_ due to the 
non-negligible contributions of the (lowest lying) even parity state. For our lattice with 
T = 40, it is sufficiently large to extract the mass of the lowest lying state for the entire 
range of 0.03 < m q a < 0.8. However, the excited state seems to suffer from the (backward 
propagating) contribution of the lowest-lying state, especially for m q a < 0.05, which yields 
relatively larger error in extracting the mass of the negative parity state. 

3. Summary and Concluding Remarks 

Our results of baryon mass spectra are summarized in Tables ^ and |2l along with the 
experimental mass spectra listed by the Particle Data Group [Ej. The empty entries in 
the last column are baryons which have not been discovered in experiments. Thus the 
baryon masses obtained in this paper can serve as theoretical predictions of lattice QCD, 
in particular, the singly- charmed, doubly-charmed, and triply-charmed baryons. Details 
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of our results including meson mass spectra will be presented elsewhere [EJ- Our results 
of baryon masses are in good agreement with experimental values, except for some of 
the negative parity baryons, in particular, A(1405) and A c (2593). This certainly warrants 
further studies on these excited baryons (e.g., with higher statistics, larger lattices, smaller 
u(d) quark masses, and incorporation of dynamical fermions). 
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Baryon mass spectra obtained in this paper 
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2237(76) MeV]. The last column is from the 


A c 


1/2+ 


2275(27) 


2285 


listings of Particle Data Group [|1], where 


A c 


1/2- 


3027(63) 


2593 


J p has not been measured for all entries. 


s c 


1/2+ 


2458(24) 


2455 


Baryon 


J p 


Mass 


Expt 


s c 


1/2" 


2882(40) 




N 


1/2+ 


958(36) 


939 


1 — <c 


1/2+ 


2478(27) 


2466 


N 


1/2- 


1603(150) 


1535 


1 — 'C 


1/2" 


2793(68) 


2790 


A 


3/2+ 


1243(74) 


1232 




1/2+ 


2677(30) 


2697 


A 


3/2" 


1764(105) 


1700 


n c 


1/2- 


3100(130) 




A 


1/2+ 


1119(38) 


1116 


s c 


3/2+ 


2522(30) 


2520 


A 


1/2- 


1854(47) 


1405 


s c 


3/2- 


3007(47) 




£ 


1/2+ 


1211(43) 


1189 


1 — 'C 


3/2+ 


2646(36) 


2645 


s 


1/2" 


1872(118) 


1750 


1 — I C 


3/2" 


2954(94) 


2815 




1/2+ 


1326(32) 


1315 


Q c 


3/2+ 


2756(32) 






1/2" 


1935(110) 


1950 


n c 


3/2" 


3224(62) 




s 


3/2+ 


1376(61) 


1385 


^cc 


1/2+ 


3522(16) 


3520 


s 


3/2" 


1931(93) 


1940 


^cc 


1/2" 


3940(47) 






3/2+ 


1539(42) 


1530 


^cc 


1/2+ 


3637(23) 






3/2- 


2039(95) 


2030 


^cc 


1/2" 


4053(57) 




n 


3/2+ 


1668(29) 


1672 


1 — 'CC 


3/2+ 


3655(20) 




n 


3/2- 


2229(50) 


2250 


1 — CC 


3/2- 


4043(22) 




^ccc 


3/2+ 


4681(28) 




^cc 


3/2+ 


3762(17) 




^ccc 


3/2- 


5066(48) 






3/2- 


4147(31) 
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